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Abstract 

The stability of photon trajectories in models of the Universe that 
have constant spatial curvature is determined by the sign of the curva- 
ture: they are exponentially unstable if the curvature is negative and 
stable if it is positive or zero. We demonstrate that random fluctu- 
ations in the curvature provide an additional stabilizing mechanism. 
This mechanism is analogous to the one responsible for stabilizing the 
stochastic Kapitsa pendulum. When the mean curvature is negative it 
is capable of stabilizing the photon trajectories; when the mean curva- 
ture is zero or positive it determines the characteristic frequency with 
which neighbouring trajectories oscillate about each other. In constant 
negative curvature models of the Universe that have compact topol- 
ogy, exponential instability implies chaos (e.g. mixing) in the photon 
dynamics. We discuss some consequences of stochastic stabilization in 
this context. 



One of the fundamental questions concerning the dynamical properties 
of photons in the Universe is whether their trajectories are stable or unsta- 
ble. This strongly influences both the images of distant objects as well as 
fluctuations in the cosmic microwave background (CMB). In cosmological 
models in which the spatial geometry of the Universe has constant curva- 
ture K the photon trajectories (geodesies) are stable if K > 0, in which 
case neighbouring trajectories oscillate about each other with a character- 
istic frequency determined by K, and exponentially unstable if K < 0, in 
which case neighbouring trajectories diverge exponentially quickly with a 
Liapunov exponent determined by \K\. 

It is obvious, however, that the Universe is not exactly homogeneous and 
isotropic: matter is not uniformly distributed but is organized into galaxies, 
clusters of galaxies, and even superclusters of galaxies [Q. Consequently, 
the spatial curvature cannot be constant, but fluctuates. Our purpose here 
is to investigate the influence of these fluctuations on the stability of the 
photon trajectories. It might be thought that random fluctuations in the 
curvature would be a source of instability. We show that this is not in 
fact the case; their influence is a stabilizing one in that the curvature is 
renormalized by an additional (positive) factor related to the fluctuation 
amplitude and length-scales. The mechanism responsible for this stochastic 
stabilization is analogous to the one that stabilizes the Kapitsa pendulum: a 
pendulum whose pivot is forced to move periodically or stochastically along 
the vertical [2j. 

Given the fact that for constant curvature models to be consistent with 
recent observations of the CMB the absolute value of the mean curvature 
must be small - the density of energy in the Universe is within about 1% of 
its critical value [H] - stochastic stabilization is capable of dominating the 
stability balance. For example, an order of magnitude estimate suggests that 
it is relevant in the present cosmological epoch. The difference between a 
constant curvature model with K close to zero and one in which stochastic 
stabilization dominates is quantifiable in terms of the stability frequency 
with which neighbouring trajectories oscillate about each other. 

Our results apply to both open and compact geometries. There has, 
however, recently been considerable attention focused on negative curvature 
models that have compact topology [U 03 H3 13 El EH E3 • Photon dynamics 
in compact spaces is recurrent and hence, by virtue of the exponential in- 
stability, strongly chaotic (e.g. mixing) when the curvature is constant and 
negative. The fact that fluctuations in the cosmic microwave background 
(CMB) have a distribution close to gaussian has been related to Berry's 
random wave model .11 for wave modes in chaotic systems [3 and 
the phenomenon of scarring in these wave modes |T2] has been linked with 
anisotropic structures in the CMB and in the distribution of galaxies [3 . 
We discuss some implications of stochastic stabilization in this context. 

We begin by deriving the appropriate form of the geodesic deviation 
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equation - the equation that governs the separation of neighbouring trajec- 
tories. Since the aim of this Letter is to establish the principle of stochastic 
stabilization for photon trajectories, the calculations we report contain the 
essential ingredients for the effect, but ignore many additional, compara- 
tively weaker phenomena present in the early Universe. With this in mind, 
we begin with an expanding isotropic homogeneous (Friedmann) cosmology 
perturbed by density fluctuations with non-relativistic velocities, neglecting 
perturbations of vector and tensor character such as gravitational waves, and 
pressure fluctuations, for example caused by relativistic neutrinos. In the 
coordinate system called the 'conformal Newtonian' or 'longitudinal' gauge 
|13j . with the spatial variables in Robertson- Walker form, such a spacetime 
is described by the metric 

ds 2 = R 2 (t){(1 + 2<D)dr 2 - ri ^/~ 2 * 2 m 2 (^ 2 + d V 2 + 

[1 + ^(x 2 + y z + z l )\ l 

(1) 

Here R(t) is the scale factor or spatial curvature radius of the Universe, 
r = J R _1 dt is conformal time, $ << 1 is the Newtonian gravitational 
potential, K is the dimensionless spatial curvature (e.g. K = — 1 corresponds 
to a hyperbolic geometry), (x,y,z) are comoving coordinates expanding at 
the same rate as the Universe, and units are chosen in which Newton's 
constant G and the speed of light c are equal to unity. The paths of photons 
through the spacetime are null geodesies described by the equation ^1] 

where n is the photon direction, and Vj_ is the gradient in comoving coordi- 
nates perpendicular to n. Note that since <3? and its derivative are small, the 
total change of n is small, and so the transverse derivative can be replaced 
by the derivative transverse to the observed (final) direction of the photon. 
The separation of two closely spaced photons propagates according to the 
geodesic deviation equation which under our approximations leads to 

d 2 a\ _ (K + 2% 2* ft ^ ((\ (3) 



dT 2 \vJ V 2 ®tv K + 2$ m J \ji 

where £ and r\ are the components of the separations orthogonal to n, and 
there are contributions arising from both the spatial curvature K and sec- 
ond derivatives of the gravitational potential in directions orthogonal to n, 
namely tidal forces. 

These tidal forces are, in principle, entirely characterized by a complete 
knowledge of the fluctuations in the matter density. We shall consider them 
to be random functions of position (and possibly time). As it moves under 
their influence, a photon will be seen to be acted on by a time-varying force 
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which fluctuates rapidly (because the speed of light is large and the length 
scale of the density fluctuations is small compared to R) and randomly 
(i.e. stochastically) with zero mean. We thus replace the spatially random 
tidal force terms in © with a time-dependent stochastic perturbation. 

In order to illustrate the qualitative behaviour of the solutions of the 
resulting class of equations, we consider first the analogous one-component 
case: 

d 2 u 

— = -{k + Af{r))u. (4) 

Here A is a control parameter and /(r) is a stochastic forcing function, which 
we take to have zero mean. This has a mechanical analogy: it describes a 
pendulum with a vertically moving pivot in the limit of small oscillations. 
In this case u denotes the angular displacement from the vertical and the 
forcing term describes the height of the pivot. This problem was studied 
in detail by Kapitsa [2]. The constant gravitational force represented by k 
(k < corresponds to an up-turned pendulum and k > to a down-turned 
one) plays the role of the smooth geometry in the case of the unperturbed 
cosmology, while the motion of the pivot corresponds to the metric pertur- 
bations induced by fluctuations in the matter density. 

We shall be interested in the case when the frequency u characterizing 
the fluctuations in /(t) is large. Stabilization has been proved in the limit 
to — ► oo using Liapunov exponent techniques (see, for example, JH])- The 
dependence of the stability on uj can be deduced by separating u asymp- 
totically into fast and slow components 0: u(t) = (u(t)) + Uf(r), where 
(. . .) denotes a local time average over scales large compared to u^ 1 . Here 
(u(t)) is the slow component and, as uj — > oo, ut is small and varies rapidly. 
Substituting into (j3J) and averaging, we have 

^ = -k(u)-A(f(r)u f (r)}. (5) 

Subtracting (J5J) from @ gives, asw^ oo, 

d 2 u f 



dr 2 



Af(r)(u(r)). (6) 



This equation can be integrated directly, treating (u(t)} as a constant, lead- 
ing to 

^ » -A(u(r)) £ f{t)dt ee -A(u(t))v(t), (7) 
where t\ is chosen so that {v ) = 0, and 

Uf{r) -A{u(t)) I v(t)dt = -A{u(t))x(t), (8) 
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Figure 1: Solutions of ® with = — 1 and /(r) given by (|T0|) . When 
j4 = 0, the solution increases exponentially, corresponding to the fact that 
the photon trajectories are unstable in a Universe with constant negative 
curvature. This is also the case when A = 20 and A = 30. When A = 
50, 60, 70 and 80 the solution is oscillatory, indicating stability. 



t 2 being chosen so that (x) = 0. Substituting into © and computing the 
integral implicit in the average by parts, we obtain 



d 2 (u) 
dr 2 



-{k + A 2 {v 2 {r))){u) 



(9) 



It follows from © that if A 2 (v 2 } < -k then u grows exponentially as 
r — ► oo, but if A 2 (v 2 ) > —k then u is linearly stabilized. This effect, and 
the accuracy of ©, are illustrated by numerical simulations, the results of 
which are represented in Figures ^ [2] and |21 In these simulations we took 
k = — 1 and 



N 



f( T ) = ^sin(w n r + 



(10) 



n=l 



with N = 100, frequencies oj n chosen at random uniformly from [120, 600] 
and phases <j) n chosen at random uniformly from (— 7r,7r]. Thus (v 2 ) = 
1 S n to^ ~ TSo an d so f° r ^ < V1440 « 38 the expectation is that u grows 
like exp(T-y/l — A 2 /1440), while for A > V1440 it oscillates with frequency 
^2/1440 - 1. 

The above analysis generalizes directly to the geodesic deviation equation 
J2J). This may be written as the stochastic matrix equation 



d 2 a 



dr 2 \rj 



K + Af n {r) Af 12 (r) \(£ 
Af 12 (r) K + Af 22 (r)J U 



11) 



4 



Figure 2: Plots of logu vs. r for the data corresponding to A=20 and A=30 
in Figure ^ The continuous black lines are best-fitting straight lines when 
r > 6. For A = 20 the gradient is about 0.85 (the theoretical value is 
0.8498. . .). For A = 30 the gradient is about 0.62 (the theoretical value is 
0.6124. . .). 




Figure 3: Plots of the fourier transform of the data for u(t) shown in Figure 
HI as a function of period 2-k/v (i.e. v is the frequency conjugate to r) 
when A = 50, 60, 70 and 80. In each case the vertical dotted line marks the 
theoretical oscillation period, 27r/y / ^4 2 /1440 — 1. 
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The variables £ and 77 can be separated into slow and fast components, as for 
u. Following the steps that led from Q to Q we find, in the high-frequency 
limit, that 



*nw - (Kl+Alv)) {W (12) 

where I is the 2x2 identity matrix and v is the matrix with elements 
Vij = f T fij(t)dt, such that (vij) = 0. The stability of the trajectories thus 
depends on the eigenvalues of the matrix i^I + ^4 2 (v 2 ); specifically, assuming 
( v n) = (^22), they are unstable if A 2 {(vn ± vu) 2 ) < —K, and if A 2 {(vn ± 
v i2) 2 } > —K they are stable. Note that v\\ and vvz are of the order of 
uj" 1 , and so the balance between the terms in the stability criterion depends 
on the size of A 2 /ui 2 \K\. In particular, since uj is considered large here, if 
K < the perturbations must be large enough to cause the curvature to be 
positive in places for stabilization to occur. This is consistent with the fact 
that the geodesies on compact surfaces where the curvature is everywhere 
negative (but not necessarily constant) are strongly chaotic. 

Let us make an order of magnitude estimate to determine the cosmo- 
logical epoch when the stochastic term is observable. Since the curvature 
is small, it will always be dominated by an observable stochastic term (see 
below). Consider that the Universe consists of a volume fraction (i?*/i?) 3 
of clusters of galaxies of size L ~ IMpc and mass M ~ 10 _4 Mpc in rel- 
ativistic units. The number of clusters within a ball of radius R is thus 
~ (R*/L) 3 . The scale factor R is set to the scale of the visible Universe at 
present, roughly 10 4 Mpc. R* is the distance scale at which voids between 
the galaxies began to form, at a redshift of order 10 (i.e. R* ~ 10 3 Mpc). 
The mean time (in terms of r) for a photon to move between clusters is of 
the order of LR 2 /B*. This gives u ~ Rl/(LR 2 ), which is large up to the 
present, but will fall to less than one in a few times the current age of the 
Universe, as the galaxies move apart. Note that uj is assumed large in the 
calculations described above. 

The tidal forces are, up to factors of 4tt, of the order of the dimen- 
sionless density MR 2 /L 3 , from Poisson's equation when a photon passes 
near a cluster; Av is typically this quantity integrated over a mean free 
time, MR 4 / L 2 R3. \K\ < 0.01 in these units. The stochastic stabiliza- 
tion will be observable if Av > 1, and will dominate the curvature if 
Av/\K\ > 1 (a weaker condition). Taking this together with the estimate for 
uj, we conclude that stochastic stabilization is potentially observable when 
(^/Mi?*) 1 / 4 < R/R # < {R^/L) 1 / 2 ; that is, for a period starting soon af- 
ter the galaxies formed and due to end a few times the present age of the 
Universe from now. The magnitude of each deflection is determined by 
M/L ~ 10~ 4 , so this effect should be observable on angular scales of arcsec- 
onds or smaller ^1]. From this perspective it should be easier to observe in 
distant quasars, which are almost pointlike, than the almost homogeneous 
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cosmic background radiation. 

The stochastic stabilization of photon trajectories has several interesting 
consequences for cosmology. In particular when K ~ it is observable in 
the characteristic frequency with which neighbouring orbits oscillate about 
each other, corresponding to a periodic focusing and defocusing of distant 
images. 

Given the recent focus on compact models of the Universe ^U] , we devote 
some concluding remarks to this case (although, as already noted, the mech- 
anism is independent of topology). The simplest models of the Universe that 
go beyond having spatially constant negative curvature incorporate static 
fluctuations. Such models are to a large extent still unrealistic, because the 
curvature is unlikely to be static over the timescales associated with pho- 
ton recurrence. Nevertheless, they illustrate most clearly the implications 
of stabilization, some of which we now list. 

As a result of stabilization, the photon dynamics will generically possess 
both regular (stable) and irregular (chaotic) components, rather than being 
fully chaotic. While wave modes in fully chaotic systems are believed to 
have a gaussian value distribution, it is well established in the context of 
quantum chaos that regular trajectories lead to a quantifiable non-gaussian 
component whose precise form depends on the size and position of the regu- 
lar regions in phase space This will in turn give rise to a non-gaussian 
component in the CMB, via the connection described in [JJE3- There is, in 
addition, likely to be a second non-gaussian component due to the photon 
trajectories close to bifurcation (i.e. close to making the transition between 
instability and stability). It was established in JH] that bifurcations give 
rise to non-gaussian fluctuations in wave modes, quantified by characteristic 
scaling exponents in the wavelength dependence of the moments of their 
value distribution in the short-wavelength limit. In situations when all of 
the generic bifurcations contribute, the moment exponents take on universal 
values, which were calculated in [EJ. These non-gaussian statistics are also 
likely to feed through to the fluctuations in the CMB in the models under 
consideration here. 

Bifurcations have a strong influence on the scarring of wave modes as 
well. In fully chaotic systems wave modes may be scarred by periodic orbits 
|12j . This effect is dramatically enhanced when the orbits undergo bifur- 
cation JH]) giving rise to superscars. It has been suggested by Levin and 
Barrow, in the context of constant negative curvature models, that scars 
may give rise to anisotropic structures in the CMB [JJ E3 • This mechanism 
would therefore be strongly enhanced by stabilization. Levin and Barrow 
[7J have also put forward the idea that the anisotropic structures asso- 
ciated with scars may show up in the distribution of matter. Regions of 
regular motion and superscars would significantly amplify this effect. In 
particular, stable islands typically have a fractal distribution, and this may 
provide an explanation, within the context of their theory, for the fractal 
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hierarchy of structures seen in galaxies and clusters of galaxies. Stable orbits 
also strongly inhibit the rate of mixing, giving rise to intermittency, and so 
may play a role in the suggested links between mixing and preinflationary 
homogenization 

We gratefully acknowledge very helpful discussions with Dr Gopal Basak 
and Professors Sir Michael Berry and Mark Birkinshaw. 
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